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Abstract 

We consider graphs made of one-dimensional wires connected at vertices and on which 
may live a scalar potential. We are interested in a scattering situation where the graph is 
connected to infinite leads. We investigate relations between the scattering matrix and the 
continuous part of the local density of states, the injectivities, emissivities and partial local 
density of states. Those latter quantities can be obtained by attaching an extra lead at the 
point of interest and by investigating the transport in the limit of zero transmission into 
the additional lead. In addition to the continuous part related to the scattering states, the 
spectrum of graphs may present a discrete part related to states that remain uncoupled to 
the external leads. The theory is illustrated with the help of a few simple examples. 

PACS : 03.65.Nk, 73.23.-b 

1 Introduction 

Thanks to the powerful experimental techniques used in mesoscopic physics during the past 20 
years, many interesting and fundamental effects have been investigated on systems that can be 
modelized by graphs : networks of wires through which an electrical current can flow. The 
most famous example is of course the Aharonov-Bohm oscillations in a coherent metallic ring 
Jl], ||, ||, |H (see also the excellent review but there have also been many other realizations 
of graphs like the recent ones devoted to the Aharonov-Bohm cage effect p, H, H|. Graphs 
which could seem at first sight to be oversimplified models for mesoscopic networks, succeeded 
in those cases to describe the interesting physical effects, which explain why they are so widely 
used in many theoretical works. 

Among all the useful concepts of mesoscopic physics the scattering approach plays a central 
role. It provides a powerful tool to study many physical quantities related to transport, noise, 
etc. A very important concept of the scattering theory is the one related to the Krein-Friedel 



relation, or Friedel sum rule [10, 11], which establishes a relation between the scattering and the 



spectral properties^. The use of this kind of relations in mesoscopic physics allows to express, 

1 In fact the idea of relating the spectral properties to the scattering properties goes back to 1937 when Beth 
and Uhlenbeck jL2| related the second virial coefficient of a gas of interacting particles to the phase shifts of the 
2-body scattering problem (see also or §77 of jwj). 
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for instance, the charge distribution in terms of the scattering properties. However it permits to 
consider only the total charge in the scattering region in an equilibrium situation. Several works 
have been devoted to build corresponding concepts to describe out-of-equilibrium situations 
and also to relate not only global but local quantities to the scattering (like the local density 
of states). This has lead to the introduction of the concepts of partial local density of states, 
injectivities, emissivities, etc [15, [l6|, 17]. Let us also mention the recent work on the relation 



between scattering and local density of states for the particular case of quasi-one-dimensional 
systems [ 



19]. 



The scattering theory on graphs has attracted the interest of many authors among which 
we can quote [20, 21, 22, 23, 24, 25, |2q| . The purpose of this article is to discuss the relation 
between scattering properties of the graph and the local quantities mentioned above. 



2 Motivation 

The local density of states (LDoS) is by definition : p(x; E) = (x \5(E — H)\ x). If we consider 
a graph Q connected to infinite leads, its spectrum is continuous. We can define the DoS of 
the scattering region (the graph) as pg(E) = fgdxp(x;E), where the integral runs over all 
the bonds of the graph (by convention, we do not include the infinite leads in what we call 
"graph G"). It is well known that this object is related to the scattering properties through the 
Krein-Friedel relation (or Friedel sum rule) . Graphs have the particularity (which does not occur 
in Id for example) that due to certain symmetries some states may remain uncoupled to the 
external leads. The presence of such localized states leads to the existence of a discrete part in 
the spectrum, superimposed on the continuous part : pg(E) = p Teg (E) + pdi s (E) where p reg (E) 
is the contribution of the stationary scattering states and Pdi s (E) = ^2 n g n ^(E — E n ) is the 
contribution of the localized states (see appendix |A|). The existence of these two contributions 



was already noticed by J. -P. Roth in a work [27] in which he extended the trace formula obtained 
for closed graphs without potential [^] to open graphs connected to semi infinite wires. What 
is rather unusual is that the continuous part and the discrete part live on the same intervals of 
energy. If the wave function is chosen continuous at the vertices, the absence of hybridization of 
the localized states with the states of the continuum occurs when some states vanish at all the 
vertices to which external leads are connected. It was shown in |2(| that the Friedel sum rule 
has to be modified since, obviously, only the continuous part is related to the scattering matrix : 
p reg (E) = IndetE + Tr {£ — £t }) where S is the scattering matrix of the graph. 

Similarly, the LDoS can be separated into a continuous part related to the stationary scat- 
tering states and a discrete part given by the localized states (see appendix |A|) : 

p(x; E) = J2 \^e\x)\ 2 + E n) K;(z)| 2 (1) 

a n j=l 

(if x ^ G, only the first term remains since the localized wave functions vanish outside the graph 
of course), tp^ (x) is the stationary scattering state corresponding to the injection of a plane 
wave at lead a. The sum over a runs over the L vertices connected to leads. The normalization 
is chosen to associate to those states a measure dE. The wave function (p n j(x) is normalized 
to unity in the graph, j being a degeneracy label. g n is the number of the localized states with 
energy E n . 

The modified Friedel sum rule is an example of a relation between a "global" quantity 
characterizing the graph (the regular part p reg (E) of the DoS pg(E)) and the scattering matrix. 
The purpose of the present work is to demonstrate some relations between local quantities such 
as the local density of states (LDoS), injectivities, emissivities and partial LDoS, and the objects 
of scattering theory (stationary states, scattering matrix) in the context of graphs. 
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In particular we will show that the first term of the LDoS is related to a functional derivative 
of the scattering matrix. The second term of the LDoS, if it exists, cannot be probed from the 
scattering properties. 



3 The scattering matrix 

We consider a graph Q being the domain of the Schrodinger operator H = — D^. + V(x), where 
Da; = d x — iA(x) is the covariant derivative. The graph is a network of B bonds joining V 
vertices (denoted by greek letters a, [3, . . .). Each bond (a/3) is identified with the interval [0, l a p] 
of R so that a scalar function ip(x) defined on the graph is characterized by its B components 
i/)( a p\(x). The Schrodinger operator acts on functions i/)(x) which are chosen, in a first step, to be 
continuous at the vertices of the graph : ifir a p\(x = 0) = ip a for all vertices j3 neighbours of a. We 
denote by tp a the value of the wave function at the vertex. Additionally we must add a constraint 
on the derivatives of ip to ensure current conservation. As soon as continuity of iJj(x) is required, 
the most general additional condition is a a pD x i/j( aj3 }(a) = \ Q ip a where the presence of the 
connectivity matrix a a/ g in the sum ensures that it runs over the neighbouring vertices of a. The 
connectivity (or adjacency) matrix describes the topology of the graph : a a p = 1 if a and (3 are 
connected by a bond, and otherwise. The notation ^(0/3) ( a ) = ^(ap){ x = 0) designates the 
value of the component at the vertex. A Q is a real parameter that affects the scattering at the 
vertex. It allows to interpolate between Neumann boundary conditions (A Q = 0) and Dirichlet 
boundary conditions (A Q = 00, which imposes ijj a = 0). We can develop an intuition of the 
role of this parameter by noting that if the vertex has coordination 2, the boundary condition 
describes a potential \ a 5(x) at the vertex. Note that the transmission amplitude through the 
vertex is 2/ (m a +\\ a /k) where m a is its coordination ; the transmission is maximized for \ a = 0. 

Among the V vertices of the graph, L are connected to infinite leads. The couplings to 
the leads can be chosen arbitrary such that we can go continuously from a situation where 
the graph is coupled to the leads to a situation where it becomes decoupled from some leads. 
This procedure introduces some discontinuity between the wave function at the extremity of 



the lead and at the vertex of the graph to which the lead is connected [25]. The scattering 
matrix S is an L x L matrix that depends on the energy E = k 2 . It can be constructed by 
manipulating matrices that encode the information about the graph (topology, potentials on the 
bonds, lengths of the bonds, magnetic fluxes, couplings to the leads) po| : 



S = -1 + 2W (M + W T W) 1 W T . (2) 
The rectangular matrix W encodes the information on the way the graph is connected to leads : 

W a/3 = w a 5 a /3 (3) 

with a G V ex t and j3 £ V, where V = {1, • • • , V} is the set of vertices and V ex t the set of vertices 
connected to leads (Card(V ex t) = L). The parameter w a S M describes the coupling between 



the graph and the lead at vertex a. Its precise definition is given in [25| : the transmission 
amplitude between the lead and the graph is 2w a /(l + w^). 

We call x a p £ [0,l a /3] the coordinate on the bond (a(3), starting from a (note that x a p + 
%(3a = la/3, the length of the bond). To describe the potentials on the bonds, we introduce two 
real functions f a /3(x a p), fp a {x a p) : the two linearly independent solutions of the Schrodinger 
equation [E + d x — Vr a p)(x)]f(x) = on the bond, satisfying boundary conditions : f a p(0) = 1, 
fapilap) = 0, //3«(0) = and fpailap) = 1- Fo r example, in the free case (V(x) = 0), we have 

f \ _ sink(l al3 -x a0 ) , f / \ _ sin kx a/3 
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The matrix M that contains all the information on the isolated graph (potential on the bond, 
topology) is : 



M aP (-E) 



Aa ^ ^ O'afi 



afi 



dx 



(a) 



afi 



+ a a /3 



df a /3 



dx 



a/3 



'a/3 



(4) 



we introduced the obvious notation f a p{ti) = / Q/ 3(a), f a p{lap) = fapiP), etc. is the magnetic 
flux along the bond. This matrix was introduced in the study of the spectral determinant of 



isolated graphs |3(J HH- Instead of encoding the information about the potential through 
the functions f a p{ x ), it can be more conveniently related to the reflection and transmission 
coefficients of each bond [25] : 



M, 



a/3 



la/3 



A a 



(1 ^*a/j) (1 ~i~ Tfia) ~\~ t a ^ t^ a 
(1 + fa.fi) (1 "I - ffia) ^a/i ^/xa 



2t 



la/3 



a 3 



(1 + + r f3a) ~ t af 3 tf3c 



These equations generalize the result known in the absence of the potential 
case we recover from (||) the well-known matrix: 



M a p = i 5 a /3 CLafi COtg kla.fi ~ CL af 3 



I a/3 

sin A;/ Q( g 



(5) 

In this latter 
(6) 



We are now ready to discuss the extraction of local information from the scattering matrix. 



4 Functional derivative of the scattering matrix 

The scattering matrix is a functional of the potential V{x) and we are now going to compute 
g^Mfi ■ As a starting point it is useful to note that if a 5 potential at x = xq is added to the 
potential, the first perturbative correction to the scattering matrix is exactly the functional 



derivative [16]: 



= nV(x) + X5(x- a*)] = E[V(x)] + X w ^\y( x )] + • • • (7) 

The advantage with graphs is that the addition of a 8 potential at x is easily implemented : 
it is done by adding a vertex of weight A at x, obtaining a graph Q x . The "weight" is the 



real parameter involved in the mixed boundary conditions introduced above [21, 32, 29]. As a 
consequence, the construction of the matrix E A of Q x in the vertex approach requires to consider 
matrices M and W of size (V + 1) x (V + 1) and L x (V + 1), respectively. We call W x the 
matrix describing the coupling of the graph to the leads when the additional vertex is at x (the 
matrix W x has only an additional column of O's compare to W) and M x the new matrix M. 
We have : 

EA = - 1+m - M> + W^ T (8) 
From (HH) we see that the matrix M x depends linearly on A : 

M X = M X + jK x (9) 

where M x is the matrix of the graph Q x=0 . This graph differs from Q only by its number of 
vertices : it possesses an additional vertex of weight A = at x. Since no scattering occurs at 
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this vertex when A = 0, we do not change the properties of the graph, but only the size of the 
matrices describing it. In the following we adopt notations such that all matrices with label x 
refer to the graph Q x = Q x=0 . The matrix K x contains only one non zero element coupling the 
vertex x to itself : 

(Kx)a/3 = S a p8 ax (10) 

where the indices run over the V vertices of the initial graph and the additional vertex at x. 
Expanding the scattering matrix S A in powers of A we get : 



The first term is the scattering matrix £ = X 



A=o 



1 1 

■ U2E W T = -1 + 2W W T (12) 

+ x M x + WjW x x M + W^W 1 ' 



and the second gives the functional derivative. Then : 



2 V g - , * r K x „ ; * < (13) 



5V(x) " k X M X + W^W X X M X + WjW x 
This expression allows to compute explicitly the functional derivative by manipulating matrices. 



5 Functional derivative and its relation to wave functions 

We now show the relation between the functional derivative of the scattering matrix and the 
stationary scattering state wave functions. The wave function at the vertices of Q x is a solution 

1 + E = W X V X (14) 
Wjil-H) = M X V X , (15) 

where ^ x is the (V+l) x L matrix gathering the wave function at vertex fi for the stationary 
scattering state ip^ (x) : the matrix element is by definition = ip^ . The index /i runs over 
the V vertices and the additional vertex x. a runs over the L vertices connected to leads. The 
wave function with the correct normalization is : 

^)( X ) = -L=4 ) («)( X ). (16) 



In the previous papers 26] we have used the notation ip^\x). Here we omit the energy 
label E to lighten the expressions. It follows from the above equations that the wave function 
at the V vertices for the L stationary scattering states is encoded in the matrix : 

** = 2 A4TW' J (17) 

For the state associated to a measure dE : 

^TSm^W™*- (18) 

These matrix expressions are useful if we want to establish a relation between the functional 
derivative of the scattering matrix and wave functions. 
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5.1 Time reversed graph 

If we consider a graph described by a matrix M({9 a p}) depending on magnetic fluxes 9 a p, it 
follows from the construction of M (see formula (48) and appendix A of p5| ) that the time 
reversed graph (with all fluxes reversed) is described by the matrix M({— 6 a p}) = M({# Q/ g}) T . 
The wave function at the vertices for the time-reversed graph is then : 



t.r. 



1 



w 1 



(19) 



Important remark : gives the wave function at the vertices for the time-reversed graph, 

which is related to the original graph by reversing all the fluxes. It should not be confused with 
the wave function describing the time-reversed motion of the electron (see also appendix |b|). 



5.2 A first result 

We can now relate the functional derivative (^) to the wave function. From the above remark 
we get : 



5T 
5V(x) 



(20) 



The matrix K x select the line in ^> x associated with x and the corresponding column in (\I'^ r ') T . 
Then the matrix elements read : 



ST. 



SV(x 



Since ipx°^ is the wave function at vertex x of the graph, we could write more elegantly : 

oV(x) 

For the scattering matrix element with indices interchanged we have obviously 

8n e 



5V(x 



^- = -2m^\x)^ a \x) 



which shows that 



r yt.r. 



5T 



/3a 



5V(x) 5V{x) 

as required by the symmetry of the scattering matrix under flux reversal. 



(21) 

(22) 

(23) 
(24) 



5.3 A compact formulation of the functional derivative 

All above formulae (13|,20) expressing the functional derivative of the scattering matrix involve 
matrices describing the graph Q x . The purpose of this paragraph is to simplify ( |l3| ) and express 
the functional derivative in terms of matrices of smaller size, related to the original graph Q. 

We suppose that x belongs to the bond (a/3) and we choose to organize the basis of V + 1 
vertices of Q x as {••• ,a,(3 \x} to help the discussion. The matrix of interest has the following 
structure : 



M X + W'W X 



f 















(M x ) aa + w 2 a 


(M x ) ax 




(Mr) 0/9 + Wff 


{M x ) Px 


\ 


■■ (M x ) xa (M x ) xf3 





A 


B 


C 





(25) 
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Since the vertex x is the only neighbour of a and (3, the part which is not written in the block 
A is precisely the same as in M + W T W. We have separated by a line the matrices into blocks 
related to the vertices of Q and the additional vertex x. The matrices W x and K x read : 



W x 



w\o) 












H 



(26) 



To express the inverse of the matrix M x + WjW x , we note that A - B{M X ) XX V C = M + W T W. 
Then : 





(M + W^W)' 1 


-(M + w T iy)- 1 J B(M x )-i 




-{M X )-^C{M + W'W)-' 





(M x + WjW x 



After a little bit of algebra, we obtain from (Fy 



2i 



W- 



1 



-K{x)K{xf- 



SV(x) k M + W T W w w M + VF T W 
where the V x V matrix couples only the vertices a and /3 : 

/ : \ 



(27) 



(28) 



K(x)K(xV 



{M X ) 2 XX 





(-^x)aa; 

V (M x ) Px j 



(Af x ) M (M x ) x /3 ) 



(29) 



We can also derive a more transparent expression of the matrix K(x)K(x)^ in terms of the 
two functions f a p{x) and fp a (x) introduced above. On the bond (a/3) the wave function is 



(30) 



where x = G [0, l a p\ measures the distance from the vertex a. The wave function of the 
time-reversed graph is obtained by changing the sign of flux (the wave function at vertices is of 



course also affected by this operation). Then we can write ^(^\(x a p) = (K{ x )^)n where 
tf(x)t = (... f a p(x a p)^ fpaMePe* ) ■ 

We have also : $^\x a p) = ((¥- r -) T K(x))^. 
If we now write 

1 <SE 



We find from ( ^2|) the new expression of the matrix involved in 

/ : \ 



(31) 



(32) 



K(x)K(x)* 







f a p(x a/3 )e 



(••• f a p(x aP )e w <™ fp a (x a p)e ie e* ) . (33) 



We can extract some nontrivial relations between the elements of M x and the functions f a p{x a p) 
and fp a (x a p) : 

{M x ) xa 



fap( x ap) 



(la: 



(Mr 



which could have been demonstrated more directly by constructing f a p{x a p) in terms of 
fxa(x ax ) and f x p(x x p). 



(34) 
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6 Partial LDoS, injectivies, emissivities and LDoS 



Partial LDoS. The partial LDoS is defined as [15, 16 1 

def 1 



4i7T 



p(a,x,P) 
It follows from (f|) that : 

p(a, x, (5) = Re 



J af3 



8 S Q/3 5 S a/3 



6V(x) SV(x) 



-•a/3 



K/3 f r(a) W^W 

Injectivies and emissivities. The injectivities are defined as 

P( X >P) = ^2p(a,x,(3) 



and the emissivities as 



They can be rewritten as 



p{a,x) =^2p(a,x,f3) . 





and 



p{x,a) 



p(a,x) 



2m \ SV(x) 



1 ( 



2ivr V^(^) 



(35) 

(36) 

(37) 
(38) 

(39) 
(40) 



5 St 



These two quantities are real thanks to the relations S+S^ g^M^ = and 

coming from the unitariry : S^S = SET = 1. 
We can now compute : 



SV(x) 

where we have used ( |122| ). Since 

mK x * x u -Efvw? 1 = i>i a) *w 



it follows that : 



1 / EtJ 5E 



2ivr V <5V(a;) 



(41) 



(42) 



(43) 



a/3 



a relation that has been demonstrated in a different context and by different means in Ref. [33|. 
From ( |39|) we finally get for the injectivities : 



p{x,a) = \^ a \x)\ 2 



(44) 



The physical meaning of the injectivities is now clear : it is the contribution of the stationary 
scattering state incoming from a to the continuous part of the LDoS. 
We proceed similarly as above : 



5V(x] 



■St = -2i7r (¥ x r -fK x V x tf = -2ivr (¥ x r -) T K x (* 



t.r.\* 
x J 



(45) 
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Then 

p{a,x) = (¥/^K x ¥^) (46) 



an 



and we finally get 

p(a,x) = \^ a \x)\ 2 . (47) 
Obviously we have recovered in the particular case of graphs the Onsager-Casimir relation 



relating injectivities and emissivities [15] : 



p(a,x;{6^}) = p(x,a;{-9^}) . (48) 

Injectivities and emissivities are related by reversing the sign of the magnetic field. 

LDoS. If we sum the injectivities or the emissivities, we get the same result, the contribution 

of the continuous part of the spectrum to the LDoS of the graph : 

E *)=£ «) = - i Tr { st ^y } (49) 



From (44) we demonstrate that : 

1 5 



2ivr 5V{x\ 



lndet£ = J]|V^(x)| 2 . (50) 



Once again, we insist that one of the interest of graphs is the fact that the functional 
derivatives involved in the injectivities, emissivities or the LDoS can be computed with algebraic 
calculations : 



■ In det E 



2ivr 6V(x) 



(51) 



= ^{^K(x)K(x)^} = -L Tr [ W _ M ^^ K^K^ — - ^ T ^ T } . (52) 
Relation between dE/dE and 8Y,/8V(x) 

If the relation (|43[) is integrated over the whole graph, we get the following exact relation (see 



[16] for the one-dimensional case and the demonstration in an appendix of [34] for graphs) : 

Graph SV(X) dE AE V ; 

If we trace this relation between matrices we obtain : 

ix WM hlietS = jE 1 ' lietE + ^ Tr { E - Ef } ■ < 54 » 



L 



Graph 



Note that in the usual formulation of the Friedel sum rule [10, 11, 14 1 only the first term 
appears. This is due to the fact that what is considered in this case is the variation of the 
total density of states (the LDoS integrated over the whole space) due to the introduction of 
a potential, whereas what we consider in the above formulae is the LDoS integrated in the 



scattering region only (the graph). This point is discussed in detail in the appendix B of [26] 
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Example : The ring 



We apply in this section the above considerations to a ring j^, ^] of perimeter I threaded by 
an Aharonov-Bohm flux 9. The upper arm (arc a) is of length l a and the lower arm (arc b) of 
length lb (we have I = l a + h)- We consider the situation without potential for simplicity. The 
ring is coupled to two leads at the two vertices 1 and 2 by arbitrary couplings described by the 
two real parameters w\ and wi (The maximum coupling is Wi = 1 ; in this case the scattering 
at the vertex is symmetric. The decoupling of the ring occurs for Wi — > 0). We call 9 a ^ the line 
integral of the vector potential due to the flux along the two arcs : 6 = 8 a + 6^. 




Figure 1: Ring penetrated by an Aharonov-Bohm flux coupled to two leads. The functional 
derivative of the scattering matrix at x determines the local quantities: partial density of states, 
injectivities, emissivities and the local density of states. 

The matrix M of the ring is : 

ie -ie a ie w b 



icot a + icotb 



M = I I (55) 

-^r~\^ lCOt a + lCOt b 

where cot„ h = cot kL h and s„ h = sin kl„ h- The matrix W is : 



From (Q) we obtain : 



where 



2 / iwf sin kl + w\ u%s a s b iw 1 w 2 (sbe lda + s a e l9b ) \ 
o \ 1W2W1 (sbe a + s a e b ) 1W2 sin kl + wfw^SaSb J 



S = s a s b det(M + Ty T P^) = 2(cos 9 - cos kl) + i(wj + w%) sin kl + w\wls a Sb (58) 



is the modified spectral determinant (the spectral determinant S{—E) is obtained by taking 
Wl,2 = 0). 

yd) . 

~~ * x9 





Figure 2: Scattering states incident from the left and the right. 
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The two components of the stationary state (1) wave function on the two arcs read 



^\x) 



'ivk S 

4 1] (x) = - 

and the stationary state (2) : 



1Wl _ c WaX/la 



lw l _ c -W h x/l h 



sin k(l — x) + e l6> sin kx — iw^Sb sin k(l a — x] 
sin k(l — x) + e l9 sin kx — iw 2 s a sin k(lf, — x) 



1U, 2 e ie a (x/i a -i) 



4 2 \x) 



'irk S 

iW2 _ c iOh(l-x/l h ) 
'irk S 



sin k(l — l a + x) + e ie sin k(l a — x) — iw\s sin kx 
sin k (I — If, + x) + e~ l6 sin k(lf, — x) — iw 2 s a sin kx 



(59) 
(60) 

(61) 
(62) 



From the choice of orientations of the arcs a and b (see figure) we see that for the two components, 
the coordinate x measures the distance from the vertex 1 (this is why we consider the component 

Mx) instead of fo(x)). Then ^ a) (0) = ^ a) (0) and %b a a \l a ) = ^(h)- 

Now we show how the local FSR is applied. The matrix M x describes the same graph with 
an additional vertex x on arc a : 



M, 



( \cot lx + \cot b -j-e i9b -ife - ^ 
icot 2x + icot b -^-e l62x 



\ 



— e" 



V 



Six 



S2x 



(63) 



icotx x + \cot 2x j 



in the basis of vertices {1,2, x}. The notations lx and x2 designate the two arcs replacing the 
arc a. In particular : I = l\ x + fax + h and = 9 b + 2x + X \- The coupling matrix is : 



and 



W x 



K x 



wi 
w 2 



1 



(64) 



(65) 



After a little of algebra, we can check that (0) gives fl2|). The wave function of the time-reversed 
graph is obtained by changing the signs of all fluxes : ^^(x) = if> {a) (x)\ ^_ g .. 

A more direct derivation follows from (|28|) . We can use more efficiently (|28|). x belongs to 
the arc a, then we have : 



K(x)K(x) 



1 / sin k(l a — x) 
e l9a sin kx 



sink(l a — x) e l9a sin kx ) 



(66) 



Using (28) we can check easily that 

1 5E 



^■ {2 \x) 



^a\x) ^>(X) 



(2), 



2ivr SV(x) 

Together with E this permits to calculate the density of states of interest. 



(67) 
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Wj ^ ^ w 2 
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Figure 3: Left figure : Ahronov-Bohm ring probed by an STM. The dashed region is forbidden 
to the electron gas at the interface of 2 semiconductors. The couplings to the contacts can be 
adjusted with the help of gate voltages. Right figure : the graph that models this situation. 



A ring out-of-equilibrium studied by STM 

We now give a physical application : a ring probed by scanning tunneling microscopy (STM). 

It is well known that if the system is at equilibrium, the tunneling current is related to the 
LDoS through the Bardeen formula. If now the device under consideration is put in an out-of- 
equilibrium situation by taking different potentials at the various contacts, which induces current 
flow, the Bardeen formula has to be generalized. We have to consider the conductances between 



the contact a and the STM tip. According to Refs. [17, 35 1 the corresponding transmission 
probabilities are 

r t i p ,a = 4:ir 2 ptip\t t ip\ 2 p(x,a) . (68) 



Here pn p is the DoS in the tip, ttip the transmission amplitude between the system and the tip 
and p(x, a) is the injectivity from contact a into the point x. 

In the case of the ring threaded by a flux 6, the injectivities were calculated above. Trans- 
mission from the tip into contact a are related to the emissivities : 



T a> ti P = 4ir 2 p tip \t tip \ 2 p(a,x) 



(69) 



whereas the modification of the conductances (transmission probabilities) of the system due to 
the presence of the tip involves the partial LDoS : 



rptip 

1 a/3 



T a/3 - 47r 2 ptip|£tip| 2 /o(a, x, [3) 



(70) 



where T a p is the transmission amplitude from lead j3 to a in the absence of the tip. 

We consider in detail the case where the ring is also weakly coupled to the two leads : 
w\,2 <C 1. Transmission sepctroscopy of a weakly coupled ring was already investigated in 
Ref. P]. In this case the scattering matrix presents sharp resonant (Breit-Wigner or Fano) 



structures. It is shown in [34| that, if the spectrum is non degenerate, the modulus of the 
stationary scattering state near an energy E n of the isolated graph is 



p{x,a) 



l#(*)f 



1 



E~E n vr (E - E r 



(71) 



where <p n (x) is the eigenstate of the isolated graph. T n ^ a = EnW^ipn^a)^ is the contribution 
of the coupling to the lead a to the resonance width T n = ^ Q T n ^ a . From the appendix || we 
also see that the emissivity is equal to the injectivity in the case of the ring p(x,a) = p(a,x), 
then the conductances are equal : T t i PjQ = T^tip- [Note that in general these two transmissions 
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are not equal but T tip Cf (#) = T at i p (— 0)]. Using the equations of appendix || and the definition 
of the PLDoS we show 

, ffl i s a/ 3T na [(E - E n ) 2 - r 2 ] + 2r n r nQ r n/ 3 2 

p(a,x,p) ~ ^ . 2 2 2 . (72) 

For the perfect ring (no potential), since the eigenstates of the isolated graph are <p n (x) = 

j_ e 2nmx/i w j^h n £ % none of the conductances depend on the position x of the tip and they 
vi 

only present a resonant structure as a function of the Fermi energy Ep : 

Ttip ' Q K £ {E F -E°)2 + ri (73) 

with E n = ( 2n7 \~ 8 ) 2 and T n ^ a = j\fE^w 2 a . As a function of the magnetic flux 0, the level E n (0) 
shifts and the conductance presents also a resonant structure as a function of the flux. If we 
increase the couplings w% 2 to the leads, the resonance peaks are broadened and oscillations in 
the injectivities are generated, and thus an x-dependence of T t i PiQ , since the original eigenstates 
of the isolated ring become strongly perturbed by the couplings to the leads. 



7 Functional derivative of X in the arc language 

We have demonstrated above how the functional derivative of the scattering matrix with respect 
to the potential is calculated. We have adopted a vertex point of view: all the matrices considered 
were matrices coupling the vertices of the graph. The vertex language is rather efficient in the 
sense that it leads to the consideration of compact matrices of the smallest possible size. However 
this approach supposes that one can introduce vertex variables, which can be achieved only if 
the wave function is continuous at the vertices inside the graph (the introduction of tunable 
couplings to the leads implies that the wave function at the end of the lead is different from 
the wave function at the vertex, however the wave function is still continuous inside the graph). 
The continuity of the wave function at the vertices implies a particular choice for the scattering 
at the vertices: the transmission amplitudes between all the leads issuing from the same vertex 
are equal, a description that may not be absolutely satisfactory in all cases. To describe the 
most general situation one has to abandon the constraint of continuity of the wave function 
at the vertices. In this case it is not anymore possible to define vertex variables and one has 
to introduce arc variables. We recall here some notions presented in |2E|. On each arc i we 
introduce an amplitude A4 arriving at the vertex from which i issues and an amplitude Bi 
departing from it (see figure). That is to say that the wave function ipi( x ) on the bond is 
matched with Aie~ lkx + B{e lkx at the extremity of the arc. 

<s ■ -— _ 

1 T\ 



Figure 4: The two amplitudes associated to the arc i. 

It is clear that we have to introduce L such couples of amplitudes, one for each external lead. 
These external amplitudes are gathered in L-column vectors A ext and £> ext . By definition the 
scattering matrix relates these amplitudes : B cxt = T,A ext . On the other hand we must introduce 
two couples of amplitudes A{, Bi per bond of the graph, i.e. one couple per arc. We gather 
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these 2B amplitudes into the column vectors A mt and B mt . Finally we group all amplitudes, 
internal and external, in two 2B + L column vectors A and B. 

The scattering by the bonds is described by a matrix R coupling reversed internal arcs : 
A mt = RB mt . The scattering at the vertices is described by a matrix Q coupling arcs issuing 
from the same vertex : B = QA. If the basis of arcs is organized as {internal arcs, external arcs}, 
the matrix Q is separated into blocks: 



Q 



Q' mt 




Q 


Q ext 



(74) 



For simplicity we suppose here that the vertex scattering matrix is symmetric, that is the 
scattering at the vertices is not influenced by the presence of a magnetic field. However it is 
straightforward to extend the results to the more general situation: one has to make a distinc- 



tion between the two off-diagonal blocks of Q. The scattering matrix reads [25]: 



S = Q ext + Q (i?t - Q^)- 1 Q T 



(75) 



We now follow the same methodology as in the vertex approach to derive an expression of 
the functional derivative of E involving arc matrices. We consider a new graph which is 
similar to the original graph Q, apart from the fact that it possesses an additional vertex at x 
of weight A (we recall that Q and Q x = Q x=0 possess the same properties). We call R x the new 
bond scattering matrix. It describes the same physics as R. The vertex scattering matrix is Q x . 
Compared to Q, Q x has an additional 2x2 block describing the scattering at the vertex x. This 
block couples the two arcs issuing from x : 



2 + ik/X V 1 1 



1 1 



1 



(76) 



We introduce the notation Q x = Q x °. This matrix describes the same physical situation as Q. 
If we expand Q x in powers of A we get for the internal part, at first order : 



(Q x i mt = QT - ^ + 



lint 



iA 



(77) 



{k x couples the internal arcs). The only non vanishing elements are those associated to the two 
arcs issuing from x: 

I 



\ 



... o 








••• 


... o 


1 


1 


••• 


... o 


1 


1 


••• 


... o 








••• 



I 





( 




\ 















1 








1 















\ 




/ 



1 10 



(78) 



The functional derivative of X is equal to the first order term of X A . We get 

6E 



5V(x) 



QT) 1k x{R x 



(79) 



This expression allows to compute the functional derivative of E in terms of matrices describing 
the graph Q x . The component of the stationary scattering state on the bond (a/3) to which x 
belongs is given by 



\^Kk^[ a(3) {x) 



1 1 • • • ) (4 - Q 



Qx 



(80) 
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With these results we can recover the relation (22) within the arc approach. 

We are now going to simplify this expression in the sense that we will relate the functional 
derivative to the matrices of the original graph Q, of smaller size than those of Q x . 

We suppose that x belongs to the bond (a(3) of Q. We show the structure of the matrices 
describing Q in a basis {other internal arcs, a(3, (3a | [external arcs}. The vertex scattering matrix 
is given by 

q = ( Q' nt II Q T 



V Q \\ Q cxt J 

where we separate by a double line the internal and external vertices, 
matrix is 

(-Roia 
q ( r a/3 t/3a 
V *a/3 r f3a 

where i? ia is the block coupling all internal arcs apart from a(3 and (3a. 

We now examine the structure of the matrices describing the graph Q x 
is replaced by two bonds (ax) and (x(3) in Q x . 



(81) 

The bond scattering 
(82) 

The bond (a(3) of Q 




Figure 5: Auxiliary bonds a = ax and b = (3x on the a/3. 



The two arcs will be also denoted by : a = ax and b = x(3 (see Fig. ||). For simplicity we 
chose to organize the basis of arcs as : {other internal arcs, a,b\a,b ||external arcs}. The vertex 
scattering matrix reads : 



Q, 



Q int 




Q T 
















Q 





Q ext 



\ 



(83) 



where the blocks of Q appear. We separate by a single line the part reminiscent of Q and the 
one added by considering Q x . The bond scattering matrix is 



Rx — 


( RoVA 






( r a \ 

V o n ) 




( t- a \ 

V o t b J 




( 


( t a o \ 
V o t- b ) 


( r- a \ 

V o n ) 



\ 



Note that the transmission and reflexion coefficients of the arcs a 
to the ones of the arc a(3 through the following relation : 



(84) 



ax and b = x(3 are related 



r a /3 tp a 
taf3 



r* 
r? 



t* 







t- 



-1 





f* 

1 1, 



(85) 
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It implies the expected relations 



t~~zt i etc. (A remark about 

-L f a f h y 



r o/3 r a + l-r'a,r h ' I.- '•,;/■/, " 

notations : t a R is the transmission amplitude from a to /3 since it is the transmission of the arc 
a/3 connecting a to /?.) To compute the functional derivative of the scattering matrix, we need 
to find the inverse of 



4-Q 



int 

x 



\ 



f* 





o tt 



\ 



■t -i 

-1 r* h 



A 


B ' 


C 





(86) 



One can see that 



A - BD~ X C = R) - Q int 
involves the matrices of Q. We have : 

Then after a little bit of algebra, (|79| ) gives : 

5H 



(B) - Q' mt )- 1 


-(B) -Q^y^BD' 1 


-D~ L C(B) - Q m )~ l 





(87) 



(88) 



SV(x) 



^Q(BJ - Q mt r l C{x){B) - Q int )" 1 Q T 



where we have introduced the matrix C(x), that couples only the two arcs a(3 on which is put 
the vertex x : 



C(x) 



( : 




\ 



X0 



1 r xa r x f} 1 r ia r x 



(90) 



We have now achieved our program since (q9|) involves matrices of the original graph Q, and not 
bigger matrices of the graph Q x (describing however the same situation as Q) with an additional 
vertex at x. 



Note that in the free case (no potential) the transmissions are t^ u 
matrix takes the simple form : 



-i9 u 



and the above 



C(x) 



( 



V 











g { 2ikx a p 
Q-ikl a f3+i6 a p 





-ikl a , 3 -i8 a 



\ 



/ 



(91) 



Following what has been done previously we try to obtain a more transparent expression 
of C(x). The starting point is to express the wave function on a given bond (a(5) in terms of 
the arc amplitudes. In the arc language, the appropriate basis of solutions of the Schrodinger 
equation on the bond [E + d^ — V^ a ^(x)]f(x) = is not anymore the functions f a (3(x) and 
fpa(x), but the couple of stationary scattering states 4> a p(x) and (f>pa{ x ) associated with the 
potential V^ a ^(x) on the bond (a(5). The function 4> a p{x) is the scattering state incoming on 
the bond from the vertex a and is matched out of the bond like : 4> a /3( x ) 
x < and (frapix) = t a/3 e ik( - x ~ la ^ for x > l a p W 



Q ikx _|_ \kx £ or 
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Figure 6: The two couples of amplitudes A a p, B a p and Ap a , Bp a related to the arcs a(3 and 
(3a. 

Then the component of the wave function on the bond reads : 

4>{aP){x a p) = B al3 <j) a /3(x a p) + Bp a <f>p a (x a p) ■ (92) 

If we introduce the 25-vector 

/C(x) t = (--- <p a fs{x a/3 ) <j)p a (x a p) ) (93) 

we can write 

tyaftM = /C(x)t B int . (94) 

We have now to specify the value of the amplitudes B a p and Bp a for the stationary scattering 
state. If we consider the scattering state ^ tl \ it is described by external incoming amplitudes 
A^ cxt whose components Ai^ cxt = -^=(5^ describe a plane wave entering the graph from the 
lead connected to the vertex fi. From the two equations A mt = RB mt and B = QA we get 

B^ int = (1 - Q int R)- 1 Q T A^ cxt (95) 

then 

$3nM = ~ Q int ^) _1 Q T ^ (M)ext = -7== k(x)t(l - Q^R)- 1 ^] . (96) 

Therefore, using (22) we obtain : 



' 1 -Q{1- RQ'^^RlCix^ix^il-Q^R)- 1 ^ (97) 



2ivr 5V(x) A-rrk 
which shows that 

C{x)R = /C(x)/C(x) t . (98) 

We have derived here non trivial relations between the reflection and transmission coefficients 
for the bonds (ax) and (/?x) on one hand, and the functions </> Q /3(x) and (j>p a {x) on the other 
hand. 

7.1 Injectivities and emissivities 

To compute the objects involved in the injectivities and emissivities, one needs to use the 
following relations : 

£ f Q - Q int ) -1 = Q* (R- Q int t)- 1 ^ (99) 

and 

(i? f - Q^y 1 Q T T) = R(R- Q int t)" 1 Qt . (ioo) 

The demonstration of these two relations uses (|75| ) and relations between the blocks of Q coming 
from its unitarity. The computation of injectivities requires an expression for : 

- ^mcj = (1 " rf^r^ww^ 1 - q^ Ar )~ i q' 1 (id) 
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and the emissivities require 



1 



2i?r 5V(x) 4irk 
7.2 Example : the ring with one lead 



Q(W - Q mt y ^(x)IC(xy (R-Q 



intt r lgt 



(102) 



We consider the ring with only one lead. This geometry has been used to illustrate the deco- 
herence introduced into a closed graph (here a ring) by the coupling to an external lead . We 
now apply to this geometry the arc approach. 




Figure 7: 



We have : 



and 



The matrix C(x) is : 



Q 



2 + w' 



1 


1 


w \ 






1 


w 




w 


w 


w' 2 J 





R 



ikl-ie 



C(x) 



e 

gifci+i# Q 



g— 2ikx Q—ikl—i9 
-ikl+i9 p -2ik(l-x) 



We have det(I& - Q 



2+ w ' 2 e M w here the modified spectral determinant is 



S = 2(cos 9 — cos kl) + \w 2 sin kl . 
After a little bit of algebra we get, using (p9|) : 
5 E 2iw 2 

t—. . = — — [sin 2 kl + 2 sin kx sin k(l — x) (cos 6* — cos M)l . 
SV{x) kS 2 1 J 

We can check that this expression coincides with the one proven above 



5V(x) 



2i?r V (x) *l>(x) 



(103) 

(104) 
(105) 

(106) 

(107) 

(108) 



by using the wave function ip(x) that has been computed before (it is given by the wave function 
of the ring with two leads ipa (x) in which we set W2 = 0) : ip(x) = ^^~ e l6x / l (sm k(l — x) + 
e~ l6> sin kx). 
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8 Conclusion 



In this article, we have discussed the local Friedel sum rule which allows to extract local infor- 
mation (the local density of states, injectivities, emissivities, and partial density of states) from 
the scattering matrix. For this purpose one has to be able to compute functional derivatives 
of the scattering matrix, which can be achieved for graphs by algebraic calculations. We have 
presented a discussion both within the vertex language and the arc language. The former is 
appropriate when we consider wave functions which are continuous at the vertices of the graph, 
whereas the latter is needed for the most general graphs with arbitrary vertex scattering. We 
have emphasized that the scattering can only give information about the states of the graph 
coupled to the external leads. If the graph possesses some states that remain uncoupled to the 
leads, this part of the spectrum (discrete part) is not probed by elastic single particle scattering. 
Local information on density of states is important for the solution of physical problems which 
involve a changing charge distribution, like non-linear transport and frequency dependent trans- 
port. Here we have given several examples to demonstrate the applications of our formalism 
and in particular have considered the experimental situation in which the local properties of a 
system are probed by scanning tunneling microscopy. 

The local density of states and its generalizations investigated here are not the only quantities 
related to functional derivatives of the scattering matrix. Off-diagonal elements of the functional 
derivatives of the scattering matrix are related to (dynamic) charge fluctuation s |36[| w hich play 
an important role in the description of dephasing in mesoscopic conductors |37|, |3£j| . Further- 
more low frequency parametric pumping of electrical conductors [39, fl(| and the fluctuations 



associated with it [|41[] can be related to functional derivatives. Thus the work presented here 
can be expected to be useful in a wide range of physical problems. 

Acknowledgement 

This work was supported by the Swiss National Science Foundation. 

A Spectrum of open graphs - Continuous and discrete parts 

We have already mentioned in the introduction that certain graphs may possess some localized 
states that are not coupled to the leads. Since such states are not manifesting themselves in the 
scattering properties, the usual state counting method from the scattering (Friedel sum rule) 
fails [26|. The purpose of the appendix is to discuss the structure of the LDoS when such a 



situation occurs. 

For each arc i of a graph, we introduce an amplitude A4 arriving at the vertex from which 
i issues, and an outgoing amplitude B{ (figure [|). Those amplitudes are related by vertex 
scattering B = Q A. The amplitudes associated with internal arcs are also related through the 
bond scattering matrix A mt = RB mt . If we eliminate B mt we get 

gT^cxtg = _ gint) A mt ( 1Q9 ) 

B cxt = Q A int + Q ext A cxt . (110) 

Two cases occur : 

(i) In general det(i?t _ Q^) ^ 0. Then at any energy k 2 , the solution of the above equations 
has components on all bonds of the graph and on the leads. All the solutions of the Schrodinger 
equation are the stationary scattering states. In this case, the FSR gives the correct information 
on the number of states of the grapfiQ. 

2 Note that det(B f - Q int ) ^ for the one-dimensional case. The Id case may be viewed as a single bond 
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(ii) For certain graphs det(i?t — Q mt ) = possesses a discrete set of solutions E = Ei,E 2 , ■ ■ ■■ At 
those energies one can construct solutions of the Schrodinger equation such that A ext = B ext = 
while the internal amplitudes satisfy (R) — Q mt )A mt = and QA mt = 0. These relations describe 
a wave function localized inside the graph. In this case, the continuous spectrum related with 
the stationary scattering states coexists with a discrete spectrum of localized states. The LDoS 
takes the form : 

9n 

p(x;E) = ^2\^\x)\ 2 +^2^25(E-E n )\^(x)\ 2 (111) 

a n 7=1 
' . ' » . ' 

continuous spectrum discrete spectrum 

for x in the graph (if x belongs to a lead, the second part of course vanishes). <p n ,j(&) is a wave 
function localized in the graph and normalized to unity and j a degeneracy label. The number 
of uncoupled states at energy E n is g n = dimKer(i?T — Q mt )\ E=En . 

We emphasize that what is unusual here is that the discrete and continous parts of the 
spectrum coexist at the same energies, due to the absence of hybridization of the localized 
states with the states of the continuous spectrum. 

The second situation may occur if the spectrum of the isolated graph is degenerate. In the 
space of the parameters of the graph (lengths, fluxes,...) this occur in a volume of measure 
zero. This "violation" of the FSR occurs for discrete values of the parameters and signals a 
discontinuous behaviour of the scattering matrix as a function of those parameters. 



A.l Example 



Let us consider the case of the ring coupled to one lead. This example has already been studied 
in [26] in the vertex language. Here we adopt an arc language needed to construct all the 
solutions of the Schrodinger equation in the case (ii) and compute the LDoS. The matrices R 
and Q are given above (103,104). 

fit _ Q"'t = | 



2+w 2 
-ikl+iO 



-ikl-W 



2 

2+w 2 



2+w 2 



2+w 2 



whose determinant is 

det(^ - Q int ) 



2 + w 2 



e~ ikl [2(cos 9 - cos kl) + iw 2 sin kl] 



(113) 



(i) If 7^ 0, the equation det(R) — Q mt ) = doesn't have a solution. The only solution of the 
Schrodinger equation is the stationary scattering state. Then the LDoS is 



p(x; E) = ip(x 



w 



sin k(l — x) + e ie sinkx 



r^xk 2(cos 9 — cos kl) + iw 2 sin kl 



(114) 



In the limit w — > 0, the wave function ip(x) gives the wave functions of the isolated ring. 
(ii) 6 = : However if the flux is zero (or a multiple of the flux quantum), 



det( J R t - Q 



int\ 



2 + w 2 



e~ ikl [2sin(A;Z/2) + iw 2 cos(kl/2)] sm(kl/2) 



(115) 



(two arcs a and o) connected at its extremities to two leads (arcs 1 and 2). The corresponding 4x4 vertex 
scattering matrix Q is decomposed into the four 2x2 blocks which can always be chosen as : Q ext = Q mt = 
and Q = Q T = 1. Then | det(J?' — Q lnt )| = 1. The FSR always counts correctly the states in Id, as it should. 
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vanishes for all energies of the isolated ring : k n = 2nir/l for n 6 N. At k ^ k n the solutions 
of ( 109|J110D give the stationary scattering state. At k = k n , ( |109 , 110| ) possesses a solution 
A mt = (1, 1) x A ext /w coupled to the lead (the scattering state), and a solution A mt = (1,-1) 
with A ext = B ext = localized in the ring. The LDoS is then : 



p(x-E) 



The first term is 



cos k(x — 1/2) 



/7rfc2sin(fci/2) + lw 2 cos(kl/2) 
InS. 



oo 

E 

n=l 



+ yj(E-k 2 r , 



- sin k n x 



[116) 



2m 6V(x) 
Continuity of £ 

To end the section, let us discuss the question of continuity of S. The scattering matrix should 
be a continuous function of the energy obviously, however it has no reason to be continuous as 
a function of parameters such as fluxes, lengths,... In the case of the ring we have : 



At zero flux 



S(A: 2 ,0) 



2(cos kl — cos 9) + iw 2 sin kl 
-2(cos kl — cos 9) + \w 2 sin kl 



-2sin(fc//2) + \w 2 cos(kl/2) 



(117) 



(118) 



2sin(/c//2) +iw 2 cos(kl/2) ' 

When we study the system through its scattering properties, there is no reason to introduce 
some arbitrary jumps of 5f(k 2 ), and the only natural choice is to impose the continuity of 5f(k 2 ) 
as a function of energy, by convention. Then, in the absence of localized states when the FSR 
holds, it is related to the integrated DoS (IDoS) M{E) of the graph by M{E) ~ ^5f(E) up to 
some oscillatory term inessential at the level of the Weyl term of the IDoS. 

In the case of the ring we see that S is discontinuous as a function of 9. For example at the 
energies kl = 2mr : 

-1 (119) 



lim lim E(k ,9) 

9^0kl->2n-K 



whereas 



lim lim£(r,0) = +1 



(120) 



The discontinuous behaviour is even more stricking on the Friedel phase which is plotted for 
different values of 9. If 9 ^ the Weyl term of the Friedel phase is 5f ~ 2kl, but if 9 = it only 
grows like 5f ~ kl [^6[ as illustrated on the figure. 



B Time reversal symmetry and some useful relations 



Note that the unitarity of (|12|) is demonstrated by using the fact that the matrix M is anti- 
hermitian p5| = — M and by using the obvious relation : 



+ 



1 



M + W T W -M + W T W M + W T W 



2W T W- 



-M + W T W ' 



(121) 



With the help of this relation we can easily find the relation between the wave function and 
the one of the time-reversed graph. From (018,1^) : 



or 



S t($*-r.)T = 



(122) 

(123) 
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Figure 8: Friedel phase 5f(k 2 ). The coupling is w = 0.25 and the length I = 1. The four curves 
corresponds to 9 = 1.5 (squares), 9 = 0.5 (triangles), # = 0.1 (diamonds) and 9 = (circles). 
The two lines are kl and 2kl. If 9 ^ there are two jumps of 2-7T per interval [2n7r, 2(n + l)7r[ at 
the energies of the two eigenstates of the isolated ring. If 9 = there is only one jump of 2tt at 
the degenerate energy. In the inset we clearly see the discontinuity of 5f(k 2 = 0) which equals 
vr for 9 ^ and for 9 = 0. 



These two last relations express how time reversal symmetry acts on the scattering states. 
Let us now discuss this point in more detail. If we consider a close graph Q characterized by 
a matrix M{^y), the graph Q tx ' obtained by reversing the magnetic fluxes is described by the 
matrix M(-y) T [see eq. (|],|5|)]. If the original graph Q has a spectrum {Eq, E±, ■ ■ ■}, the time 
reversed graph Q l ' r ' has the same set of energies : 

E n ({9^}) = E n {{-9^}) . (124) 

This very general feature is easily demonstrated for graphs since the spectrum is given by 
det M{—E n ) = 0. Note that stricktly speaking the energies are not in general even functions of 
the magnetic field. It is only the full spectrum which is invariant under time reversal symmetry, 



however we do not change the labels on both sides of equation (124) to simplify the discussion. 

If (p n {x) is an eigenstate of Q, then the basic fact that the corresponding eigenstate of Q t - r - 
is its complex conjugate, 

<p t n r -(x) = <p n (x)*, (125) 

can be also easily recovered. Indeed, if we call (p n the F-column vector gathering the values of 
the wave function at the nodes, <p£ = (ip n ,i> • • • > fn.y), it is a solution of M(—E n )ip n = 0. Since 
M(-E n ^ = —M(—E n ) we see that the eigenvector of M(-E n ) T is ip* n . 

Now we examine the connected graph and its stationary scattering states. What is the 
relation between the scattering states of Q and Q t,r ' ? Equation ( |122j ) relates the values of tp 
and tp t ' r ' at all the vertices. Since we can always introduce an additional vertex on a bond, it 
can be written as : 

j2K^ Lr - {a) (x) = ^Hxr . (126) 
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If the graph is weakly coupled, the scattering matrix has a resonant structure : 



~ - da/3 + 



E — E n + \T n 

2|,„ („,\\2 



(127) 



where T n = T n;0l with T n ^ a = y/E^,w^\ip n (a)\ 34]. The positions of the resonances and 
their widths remain unchanged by reversing the fluxes. The scattering states behave like [34] : 



1 iEn w a (p n (a)* 



(128) 



E~E n \/tt E - E n + \T n 
near a resonance E n , and the scattering states of Q l ' r ' behave like : 

It.r.(a), \ 1 i-En /4: w a ip n (a) 

E~E n \/TT E - E n + ll n 

From the last equations we see that we must take care of the fact that ^) l ' r ' ^ a \x) / (x)* . 



(129) 
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